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Abstract. We analyzed the spectra and light curve 
of the peculiar type II-P supernova 1997D to recover 
ejecta parameters. The optimal hydrodynamical model of 
SN 1997D, which meets observational constraints at the 
photospheric epoch, suggests a low explosion energy of 
about 10^*^ erg, ejecta mass around 6Mq, and presuper- 
nova radius near 85Rq. We confirm the previous result 
by Turatto et al. (199q) that the ejecta contain a very 
low amount of radioactive ^'^Ni (« O.OO2M0). Modelling 
the nebular spectrum supports the hydrodynamical model 
and permits us to estimate the mass of freshly synthe- 
sized oxygen (0.02-0.07 Mq). Combined with the basic 
results of stellar evolution theory the obtained parameters 
of SN 1997D imply that the progenitor was a star from the 
8-12 Mq mass range at the main sequence. The fact that 
at least some progenitors from this mass range give rise to 
core-collapse supernovae with a low kinetic energy (« 10'^*' 



K*" ■ erg) and low amount of radioactive ^^Ni (« 0.002^/©) has 

Jv>( , no precedent and imposes important constraints on the 

^ ■ explosion mechanism. We speculate that the galactic su- 

5^ ; pernovae 1054 and 1181 could be attributed to SN 1997D- 

like events. 
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1. Introduction 

The Type II supernova (SN) 1997D discovered on Jan. 
14.15 UT (De Mello & Benetti |l997| ) is a unique event 
characterized by extremely low expansion velocity, low lu- 
minosity, and very low amount (« 0.002Mq) of radioac- 
tive ^^Ni (Turatto et al. 1998). An analysis of the observa- 
tional data led Turatto et al. (1998) to conclude that they 
caught SN 1997D around day 50 after it had exploded as 
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a red supergiant with the mass of 26 Mq and radius of 
Rq « 3OOi?0. The derived ejecta mass is M « 24A/0 and 
kinetic energy is i? ~ 4 x 10^° erg. They propose a sce- 
nario in which the low ^^Ni mass in SN 1997D is caused 
by a fall-back of material onto the collapsed remnant of 
the explosion of a 25-40 Mq star. An exciting implication 
is that SN 1997D might be accompanied by black hole 
formation (Zampieri et al. 199S). 

Here we present arguments for an alternative view on 
the origin of SN 1997D, which in our opinion was a de- 
scendant from the low end of the mass range of core- 
collapse supernova (CCSN) progenitors. The problem we 
attempted to solve first was to find the hydrodynamical 
model, which could reproduce the light curve, the velocity 
at the photosphere, and the line profiles of major strong 
lines (Section H). We emphasize the importance of the line 
profile analysis, since it provides a robust information on 
the velocity at the photosphere. The latter is a crucial 
parameter for constraining hydrodynamical models. Un- 
expectedly for us it turned out that Rayleigh scattering in 
SN 1997D is significant and may be used as a powerful di- 
agnostic tool. The emphatic role of Rayleigh scattering in 
this case is related to the low energy-to-mass ratio [E/M] 
of SN 1997D (Turatto et al. |1998D , which results in the 
higher than normal density at the photospheric epoch for 
SN II-P. A combination of hydrodynamical modelling and 
robust analysis of spectra at the photospheric epoch per- 
mitted us to impose tight constraints on E, Af , and Rq of 
SN 1997D. 

In addition, we analyzed nebular spectra of SN 1997D 
using a nebular model (Section ^. To make such an anal- 
ysis as secure as possible we first checked the model taking 
advantage of the well studied SN 1987A at a similar epoch. 
We found modelling the nebular spectrum of SN 1997D 
beneficial in discriminating between low and high-mass 
options for the ejecta. To our knowledge the present pa- 
per is a first attempt of a SN II-P study to simultane- 
ously make use of all data: light curve, photospheric and 
nebular spectra. Some implications of the low mass and 
kinetic energy of the SN 1997D ejecta for the systematics 
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Fig. 1. Synthetic (thick hne) and observed (thin hne) 
spectrum of SN 1997D of Jan. 17 1997. All the strong lines 
are indicated. Note that the model overproduces emission 
in Na I Di,2 and Ba II 6142 A lines. 



of CCSN, explosion mechanism, and galactic population 
of supernova remnants (SNR) are discussed in the final 
section. Below we adopt for SN 1997D the dust extinc- 
tion Ab = 0.0 mag and the distance 13.43 Mpc following 
Turatto et al. 



(199S). 



2. Photospheric epoch 

2.1. The velocity at photosphere 

According to general results of hydrodynamical simula- 
tions of SNe II-P the radiation cooling of the expanding 
envelope at the plateau phase proceeds in a specific regime 



of the cooling recombination wave (Grassberg et al. 1971) 



As a result the photosphere in SN II-P resides at the well 
defined jump between the almost completely recombined 
(ionization degree < 10~*) transparent atmosphere and 
the fully ionized sub-photospheric layers of high opacity. 
The velocity at the photosphere determined from the ob- 
served scattering line profiles during photospheric epoch 
(plateau) thus gives us a position of the cooling recombi- 
nation wave and therefore is of vital importance for con- 
straining parameters of the hydrodynamical model. 

To measure the photospheric velocity in the Jan. 17 
spectrum of SN 1997D we concentrated on the 5600-6700 
A band which contains strong, clearly-cut spectral lines 
(Fig. g) of H, Na I, Ba II, Fe II, and Sc II. Most of them are 
well observed in other SNe II-P. However, due to the low 
expansion velocity it is possible to distinguish here some 



spectral features never observed before, e.g. Sc II 6605 A 
line (Fig. W. A Monte Carlo technique used for modelling 
the spectrum (Fig. |^) suggests an absorbing photosphere 
and a line scattering atmosphere (Schwarzschild-Schuster 
model). In total 19 lines are included for this spectral 
range. The Sobolev optical depth was computed assum- 
ing the analytical density distribution in the envelope 



p = Pq[1 + [v/v]^)'' 



(1) 



which corresponds to a plateau at velocities w < Wk and a 
steep slope p oc v^" (n « 8) in the outer layers at u > Wk- 
Parameters po and Wk are defined by the ejecta mass M, 
kinetic energy E^ and index n. The case shown in Fig. ^ 
is characterized by M — 6Mq, E = 1.2 x 10^" erg, and 
n = 8, although one may easily fit the spectrum using 
higher mass and higher energy. Since we do not solve the 
full problem of radiation transfer in ultraviolet we adopt 
that metals are singly ionized and find level populations 
assuming appropriate excitation temperature (4200 K). 
For the standard abundance we assume here, Sc II 6605 A 
line is too strong; therefore its abundance is reduced by a 
factor of two. It may well be that the odd behavior of Sc II 
line reflects different excitation conditions for Sc II and 
other metals rather than abundance pattern. We found 
that hydrogen excitation has to be cut beyond v = 1400 
km s^^ to prevent washing out the 6500 A peak. Close to 
the photosphere within the layer Av = 300 km s~^ the net 
emission in Ha is comparable to its scattering component. 
We simulated this emission assuming the line scattering 
albedo greater than unity. 

In spite of its simplicity the model is appropriate for 
the confident estimate of the photospheric velocity which 



was found to be Vp = 900 km s ^ (Fig. [I]) with an un 
certainty less than 100 km s^^ 



The value of Vp = 970 



km s^^ reported by Turatto et al. ( 1998] ) is consistent 
with the above estimate. Our choice was a compromise 
between two possibilities: (1) washing out many observed 
features in the spectrum if Vp > 1000 km s"^, and (2) 
producing significant excess in emission components for 
Na I Di,2 and Ba II 6142 A hues if Vp < 800 km s'^ The 
value Vp — 900 km s~^, although optimal, still leads to 
some extra emission in Na I Di_2, Ba II 6142 A and Fe II 
6249 A lines (Fig. |l|). Preliminary analysis indicated that 
this drawback of the model may be overcome, if Rayleigh 
scattering on neutral hydrogen is taken into account. 

2.2. Rayleigh scattering effects 

Rayleigh scattering on neutral hydrogen in the optical 
dominates over Thomson scattering at an extremely low 
ionization degree x < 0.001 which is the case for SN II-P 
atmospheres at the photospheric epoch. To get an idea 
of the role of Rayleigh scattering in the spectrum of 
SN 1997D we adopt the analytical density profile given 
by Eq.(|l|) with a power index n — ^. Let us first estimate 
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ks is the scattering coefficient. In the approximation of the 
isotropic scattering the plane albedo reads 




Fig. 2. Influence of Rayleigh scattering and diffuse re- 
flection from photosphere on the resonance line profile. In 
parentheses are given the thermalization parameter in the 
sub-photospheric layers and the Rayleigh optical depth in 
the atmosphere. 



the Rayleigh optical depth tr using the cross-section by 
Gavrila (1967) and assuming conditions of the atmosphere 
of a normal SN II-P (e.g. SN 1987A) for two extreme cases: 
completely mixed and unmixed envelopes. Assuming for 
SN 1987A £; = 1.1 X lO^^ erg, M = 15Mq, a helium/metal 
core mass Mr. = A.2Mq, and v-p — 2600 km s~^ at the 



age t = 50 d ( Woosley |1988| ; Shygeyama & Nomoto |1990| ; 
Utrobin 1993) one gets at the wavelength 6142 A (Ba II 
line) a Rayleigh optical depth of 0.07 (mixed) and 0.1 (un- 
mixed) . 

SN 1997D is essentially different in that respect. 



Adopting the ejecta model by Turatto et al. (1998), viz. to- 
tal mass M = 24Mq, helium/metal core mass M^ = QMq, 
v-p = 900 km s^^ at the expansion time i = 50 d (the 
epoch of Jan. 17) one finds the Rayleigh optical depth in 
the range 1.3-1.8 at A = 6142 A, more than one order 
of magnitude exceeding that in normal SNe II-P. For the 
model with parameters scaled-down by a factor of four 
(M = 6Mq, helium/metal core mass Mc = I.5M0), one 
obtains 0.33 < tr < 0.44. Our study showed that such 
values cannot be ignored in modelling line profiles. 

Moreover, to treat Rayleigh scattering in an adequate 
way one has to abandon the assumption of a fully absorb- 
ing photosphere and instead include a diffuse reflection of 
photons from the photosphere. We describe the diffuse re- 
flection by a plane albedo A{fi, e) which is a function of 
cosine fi of incident angle and thermalization parameter 
e = ki^/{kii -j- kg). Here fca is the absorption coefficient and 



-4(A*,e) = l-0(A*,e)\/e 



(2) 



where the functi on 0(/ i, e) is defined by the integral equa- 
tion (cf. Sobolev |1975D 



0(^,e) 



-(l-e)^(/)(/i,e) / ^-^dm, 

2 Jq fi + fii 



(3) 



which was solved numerically to create a table of A(fj,, e). 

In the absence of Rayleigh scattering, non-zero albedo 
for e — 0.3 slightly (by 4%) increases the intensity of the 
emission component compared to the purely absorbing 
photosphere (Fig. g). The difference obviously becomes 
larger for smaller thermalization parameter. Rayleigh 
scattering significantly decreases the emission component 
due to backscatter and subsequent absorption of photons 
by the photosphere in the case of e = 0.3 and tr = 1. 
Another effect of Rayleigh scattering is washing out of the 
absorption trough by continuum photons drifted from blue 
to red; this effect is especially pronounced for weak lines 
and is of minor importance for strong lines. This mod- 
elling shows how the emission excess in Na I and Ba II 
lines (Fig. |l|) may be suppressed. 

Apart from Rayleigh scattering and diffuse reflection 
by the photosphere we made two other essential modiflca- 
tions to our Monte Carlo model of line formation. First, we 
took electron scattering into account. The electron density 
distribution is recovered from the Ha line profile using a 
two-level plus continuum approximation. Second, we cal- 
culated the population of three lowest levels of Ba II using 
the observed flux in the spectrum on Jan. 17. This approx- 
imation is fairly good in analyzing the blue side of the ab- 
sorption trough of the Ba II 6142 A line. We adopted the 
standard barium abundance (Grevesse & Sauval 1998 ) and 
the Ba II fractional ionization n(BaII)/n(Ba) = 1. The 
latter seems to be a good approximation for the outer lay- 
ers of SN 1987A at the stage when strong Ba II lines are 
present (Mazzali et al. 1992). 

With the modified Monte Carlo model the synthetic 
spectrum is calculated for two relevant cases: a high-mass 
model with parameters M = 24M0 and Mc = 6Mq 
(Fig. m) and a low-mass model with parameters M ~ 
6Mq and Mc = 1.5M0 (Fig. |b). Note, that both mod- 
els have the same photospheric velocity Vp — 900 km s^^ 
and the same ratio E^q/M — 1/6, where £^50 is the ki- 
netic energy in units of 10^° erg and M in Mq. Com- 
plete mixing, which implies a minimum Rayleigh optical 
depth, gives tr = 1.3 and 0.33 for high and low-mass mod- 
els, respectively. The thermalization parameter e in sub- 
photospheric layers is 0.35 and 0.24 for high and low-mass 
models, respectively. In the high-mass model Rayleigh 
scattering suppresses the emission components of Na I 
Di^2, Ba II 6142 A, and Ba H/Fe II peak at 6500 A down 
to an unacceptably low level (Fig. gja). The low-mass case 
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Fig. 3. Synthetic spectra (thick Unas) with Rayleigh scat- 
tering and diffuse reflection from the photosphere taken 
into account, and observed spectrum of SN 1997D on Jan. 
17 (thin hne). The high-mass case (a) shows strong smear- 
ing of emission components, while the low-mass case (b) 
with moderate Rayleigh optical depth fits observations 
better. 



fits the observations fairly well (Fig. Ob) . Computations of 
spectra for different values of Rayleigh optical depth led 
us to conclude that the tolerated upper limit is 0.6. Yet 
the Rayleigh optical depth cannot be lower than « 0.3, 
otherwise the emission in the Na I Di_2 and Ba II 6142 
A lines becomes too strong. We find an optimal value is 
tr = 0.45 with an uncertainty of about 0.15 for the Jan. 
17 spectrum. 

2.3. Diagnostics of ejecta mass and kinetic energy 

The observational limitations upon Rayleigh optical depth 
in the atmosphere of SN 1997D may be combined with 
the restriction on the density in the outer layers imposed 
by the blue absorption edge of Ba II 6142 A in order to 
constrain ejecta mass and kinetic energy. The idea may 
be illustrated using a toy model, in which the supernova 
envelope is represented by a homogeneous sphere with the 
boundary velocity vq. Given the photospheric velocity and 
Rayleigh scattering optical depth one finds the product 
pvot, whereas the blue edge of the Ba II 6142 A absorption 
gives the outer velocity vq. For a specific phase t one gets 
then ejecta mass M = {47r/3)p(vot)^ and kinetic energy 
E = {3/10)Mvl 

In practice we used a more realistic density profile 
given by Eq.dlf) with a power index n = 8. In this case, 
likewise for the simple model considered above, one can 
find, accepting a certain ejecta mass, the corresponding 
value of the kinetic energy compatible with the blue edge 
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Fig. 4. The 'mass-energy" diagram for SN 1997D. Solid 
lines bound the " barium" strip derived from the observed 
Ba II 6142 A absorption; dashed lines are loci of the con- 
stant Rayleigh optical depth with indicated values (see 
Section 2.3). 



of the Ba II 6142 A absorption in the SN 1997D spectrum 
on Jan. 17. Again, we adopted a standard barium abun- 
dance with the Ba II ion as the dominant ionization state. 
Variation of the model mass under the condition that the 
Ba II 6142 A absorption is reproduced results in the cor- 
responding variation of the kinetic energy. Taking into ac- 
count uncertainties of the Ba II 6142 A absorption fit we 
found a region of allowed parameters ( "barium" strip) in 
the mass-kinetic energy (M — E) plane (Fig. Q). The lower 
and upper limits of Rayleigh optical depth, 0.3 and 0.6, 
respectively, produce another strip of allowed parameters 
("Rayleigh" strip) in this plane. The overlap of "barium" 
and "Rayleigh" strips gives a tetragonal region where the 
ejecta mass and kinetic energy of SN 1997D are confined. 
One sees that optimal values of ejecta mass should reside 
around M ~ 6Mq , while kinetic energy should be close to 
E ~ 10^° erg. 

The suggested diagnostics, unlikely useful for ordinary 
SNe II-P, proved efficient for constraining parameters of 
SN 1997D. A warning should be kept in mind that a cos- 
mic barium abundance was assumed here. This may in 
general not be the case since SN 1987A demonstrates that 
barium overabundance in SNe II-P may be as large as a 
factor of two relative to the cosmic value (Mazzali et al. 
1992). If barium abundance in SN 1997D ejecta is twice 
the cosmic value, then the "barium" strip in the M — E 
plane has to be shifted down by a factor « 1.3 towards 
lower values of kinetic energy. It is remarkable that this 
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Fig. 5. Chemical composition of ejecta adopted for the 
hydrodynamical model of SN 1997D. The hydrogen-rich 
envelope is mixed with the He shell in the inner 4.5 Mq 
of ejecta. Note that the ejecta mass does not include the 
collapsed core with a barionic mass of 1.4 Mq. 



Fig. 6. Calculated light curve (thick line ), observed pho- 
tometric data (open circles), and uppe r limi ts (v-like sym- 
bols). Data are from Turatto et al.( 1998 ). Dotted line 
shows the calculated V band luminosity assuming the 
spectrum energy distribution as in the observed spectrum 
at t « 150 d. 



diagnostics does not depend on the supernova distance. 
However there is a weak dependence on reddening via the 
colour temperature determined from 4500 A/6140 A flux 
ratio which affects the Ba II excitation. Unaccounted red- 
dening leads to the overestimation of the mass obtained 
from the Ba II line. 



2.4^. Light curve 

The light curve of SN II-P during the plateau phase is 
determined by the ejecta mass M, kinetic energy E, pre- 
SN radius i?o, the structure of the outer layers, the ^^Ni 
mass and its distribution, and the chemical composition of 
the envelope (Grassberg et al. |l97l| ; Utrobin |l989| , |l993D . 
The ^^Ni mass in SN 1997D is reliably measured by the 
light curve tail. The structure of outer layers of the pre-SN 
may normally be recovered from the initial phase of the 
light curve, which was unfortunately missed in the case of 
SN 1997D. Therefore we used a standard pre-SN density 
structure with the polytrope index of three, though mod- 
els with other density structure were also tried. The abun- 
dance of the deeper part of the envelope, e.g., the transi- 
tion region between the H-rich envelope and metal/helium 
core affects the final stage of photospheric regime and may 
therefore be probed by the light curve at the end of the 
plateau phase. In general, the parameters M, E, and Rq 
then may be found from the plateau phase duration, lu- 
minosity at plateau phase (e.g., in V band), and velocity 
at the photospheric level. In a situation when the plateau 
phase duration is unknown the optimal Rayleigh optical 
depth in the atmosphere (tr « 0.45 ± 0.15) provides the 



missing constraint. The description of the radiation hy- 
drodynamics code used for supernova study may be found 



elsewhere (Utrobin |1993| , [1996|). 

An extended grid of hydrodynamical models of SN 
1997D led us to the conclusion that requirements imposed 
by the V light curve, velocity at the photosphere, and 
Rayleigh optical depth are consistent with those estimated 
above from the M-E diagram. The optimal hydrodynami- 
cal model is characterized by the following parameters: the 
ejecta mass M = 6M0, kinetic energy E = 10^° erg, and 
pre-SN radius Rq = 85Rq . To prevent the emergence of a 
luminosity spike at the end of plateau phase and to explain 
the narrow peak of the Ha emission in the nebular spec- 
trum on day « 300, we suggest mixing between the helium 
layer and the H-rich envelope (Fig. ||). The adopted he- 
lium/metal core mass before mixing is M^ — 1.5Mq. With 
0.002 Mq of radioactive ^^Ni this choice of parameters re- 
sults in a V^ lig ht cur ve, which fits the observational data 
(Turatto et al. 1998) and is consistent with the observa- 
tional upper limits by Evans at early epochs (Fig. ra) . The 
velocity at the photosphere in this model is 830 km s~^ in 
agreement with that found from the spectrum synthesis. 
Remarkably we obtained the same E/M ratio as Turatto 
et al. (1998). In our model the first spectrum on Jan. 17 
corresponds to the epoch of 46 days after the explosion in 
good agreement with the 50 days found by Turatto et al. 



( |1998D . 

The diffusion approximation used in the hydrodynam- 
ical model breaks down at the transition from the plateau 
to the radioactive tail about < « 65 d. To reproduce the 
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Fig. 7. Synthetic spectrum for the hydrodynamical model 
(thick hnc) and observed spectrum on Jan. 17 (thin hnc). 



tail, we translated the bolometric luminosity computed in 
the hydrodynamical model into the V band luminosity us- 
ing the two assumptions about the spectrum of escaping 
radiation at the tail phase. The first one admits that the 
spectrum is black-body with the constant effective tem- 
perature calculated at i = 65 d. This gives somewhat 
higher V luminosity compared to observations at the tail 
stage (Fig. ||). An alternative approach assumes that the 
spectrum of escaping radiation during the tail phase is the 
same as in the observed nebular spectrum at t ~ 150 d 



(Turatto et al. 1998). The latter assumption is more real- 
istic and provides a good fit to observations (Fig. ^. This 
agreement justifies the adopted ^^Ni mass of 0.002 Mq 
originally obtained by Turatto et al. ( 1998| ). 



The envelope structure computed in the hydrodynam- 
ical model was then used to recalculate the synthetic spec- 
trum in a way similar to that described in Section 2.2 . The 
model spectrum agrees well with the observed spectrum 
on Jan. 17 (Fig. M). Of particular importance is an excel- 
lent fit for the emission component of the Na I 5889, 5896 
A doublet which is free of blending, thus being a reliable 
probe for Rayleigh optical depth in the atmosphere. 

Analyzing hydrodynamical models with different sets 
of input parameters gives us confidence that the envelope 
mass and the pre-SN radius are determined with the un- 
certainty of about 1 Mq and 10 i?©, respectively. There- 
fore, we estimate the ejected mass as 6 ± IMq with the 
invariant ratio E^q/M = 1/6 and the radius of pre-SN as 
85 ±10 Rq. 

A dust extinction in the host galaxy (NGC 1536) can- 
not be ruled out. It is unlikely, however, significant since 
the galaxy is nearly face-on. With some dust extinction (if 
any) the kinetic energy and/or the pre-SN radius should 
be increased accordingly. For instance, the dust extinction 



Ay = 0.1 mag suggests the 13% increase of the kinetic en- 
ergy 

3. Nebular phase 

3.1. Nebular model 

The high quality late time spectrum of SN 1997D at the 
nebular epoch t 



300 d (Turatto et al. |1998|) gives us an 
opportunity of a complementary test for the ejecta model. 
Our goal here is to reproduce all the strong lines observed 
in the spectrum, viz. Ha, [O I] 6300, 6364 A, and [Ca II] 
7291, 7324 A using the density distribution of the hydro- 
dynamical model. We assume that the ejected envelope 
consists of two distinctive regions: a core and an external 
H-rich envelope. The core in the nebular model is a macro- 
scopic mixture of radioactive ^^Ni, H-rich matter (compo- 
nent A), He-rich matter (component B), 0-rich matter 
(component C), and the rest of metals, e.g., C, Ne, Mg 
(component D). The latter does not contribute noticeably 
to the lines we address. The density of the A, B, and D 
components is equal to the model local density p, while 
the 0-rich matter of density po may be clumpy with the 
density contrast xo — Po/ P- Voids arising from the oxy- 
gen dumpiness are presumably filled in by the ^^Ni bubble 
material. 

The average gamma-ray intensity was calculated us- 
ing a formal solution of the transfer equation with known 
distribution of ^^Ni and assuming the absorption approx- 
imation with the absorption coefficient k = 0.03 cm^ g~^. 
The fraction of deposited energy lost by fast electrons on 
heating and ionization of hydrogen, helium, and oxygen 
was taken from Kozma & Fransson (1992). The rate of 



nonthermal excitation and ionization of helium in the H- 
rich matter was added to the hydrogen ionization rate 
to take account of hydrogen ionization by UV radiation 
produced by helium nonthermal excitation and ionization. 
Due to this process the Ha intensity is insensitive to the 
He/H ratio. The photoionization of hydrogen from the sec- 
ond level by hydrogen two-photon radiation and by the 
Balmer continuum were taken into account as well. The 
Balmer continuum radiation consists of the recombination 
hydrogen continuum and the rest of ultraviolet radiation 
created by the radiation cascade of the deposited energy 
of radioactive '"^^Co. This additional component of Balmer 
continuum radiation was specified assuming that a frac- 
tion p of the deposited energy is emitted in Balmer con- 
tinuum with the spectrum j^ cc v~^ . We adopted p = 0.2 
according to estimates by Xu et al. ( |1992[ ) for SN 1997A 
at the nebular epoch. 

In thermal balance only the principal coolants are in- 
cluded: hydrogen lines, C I 2967 A, 4621 A, 8727 A, 9849 
A, Mg II 2800 A, [O I] 6300, 6364 A, [Ca II] 3945 A, 
7300 A, Ca II 8600 A, and Fe II lines. For the sake of 
simplicity the total Fe II cooling rate of permitted and 
semi-forbidden lines is assumed to be equal to the cooling 
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Table 1. Parameters of nebular models 



Model 


M 


^50 


Z/Zq 


Vc (km s ^) 


Xo 


TM 


15 


11 


0.4 
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5.5 


Ml 


6 


1 
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600 


1.2 


M2 


6 


1 


1 


600 


1 


M3 


24 


4 


0.3 


780 


1 


Model 


Mc 


Mn 


Mhc 


Mo 


Aimct 


TM 


4.4 


2 


0.7 


1.2 


0.5 


Ml 


0.64 


0.12 


0.45 


0.035 


0.035 


M2 


0.64 


0.12 


0.45 


0.035 


0.035 


M3 


4.83 


0.63 
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Fig. 8. Computed (thick line) and observed (thin line) 
emission lines in the SN 1987A spectrum on day 339 (Pun 



et al. 1995) 



rate of one Mg II 2800 A line. However, unlike for the real 
Mg II line coUisional saturation is omitted to allow photon 
branching in Fe II lines. Cooling via the excitation of Fc II 
forbidden lines is represented by [Fe II] 8617 A and [Fe II] 
4287 A lines, which are the most efficient coolants for the 
relevant temperature and electron density. We include also 
adiabatic cooling; it is important in the outer region of the 
hydrogen envelope. Metals with a low ionization potential 
(Mg, Ca, Fe) are assumed singly ionized. 

With a specified density distribution and ^®Ni mass 
the primary fitting parameter is the velocity at the core 
boundary Vc, which affects line intensities via the mass of 
the mixed core Ale exposed to the intense gamma-rays. 
The amount of matter in components A, B, and C should 
then be determined from the spectrum fit. 



3.2. Model test 

Before applying the nebular model to SN 1997D it is in- 
structive to compute the nebular spectrum of the well 
studied SN 1987A. We used the CTIO spectrum corrected 



for redden ing at the epoch 339 days (Phillips et al. |1990 ; 
Pun et al. 1995| ). The primordial metal abundance (Z) is 



assumed to be 0.4 solar. The density distribution in the 
envelope is approximated by Eq. (yj) with n = 8.5. Com- 
promise values of ejecta mass Al — 15AIq and kinetic 
energy £■ = 1.1 x 10^^ erg are adopt ed (W oosley 1988 ; 



Shygeyama & Nomoto 1990 ; Utrobin 1993). Apart from 
the ^^Ni mass (0.075 Mq), we specify the amount of met- 
als in the core M^^t = 0.5Mq (component D), in line with 
the expectations for an 18-22 Mq progenitor (Woosley & 
Weaver 1995; Thielemann et al. 



1996) 



A satisfactory description of line profiles and intensi- 
ties of Ha, [O I] 6300, 6364 A, and [Ca II] 7291, 7324 A 
(Fig. H) is obtained with the test model (TM) for a sound 
choice of parameters (Table Q). The table gives ejecta mass 
(M), kinetic energy (i^so), the primordial-to-solar metal 
abundance ratio {Z/Zq), velocity at the outer boundary 
of the mixed core (wc), oxygen density contrast (xo)j core 
mass (Mc), and other core components, viz. H-rich matter 
(Mh), He-rich matter (Mhc), 0-rich matter (Mq), and 
metals (Mmot)- AH masses in Table |l] are given in solar 
masses. The amount of H-rich and He-rich matter in the 
mixed core inside 2000 km s^"'^ {"^Mq and O.7M0, respec- 
tively) are in good agreement with values advocated by 
Kozma & Fransson (1998). The rest of newly synthesized 
helium (« 1-Mq) is presumably mixed with the H-rich 
component. The oxygen mass and density contrast were 
found from best "eye-fit" of flux of the [O I] doublet. The 
value Xo = 5-5 corresponds to the oxygen filling factor 
K, 0.045, a value earlier found by Andronova (1992). 

The oxygen mass estimate is hampered somewhat by 
the uncertainty arising from the poorly known fraction 
of oxygen cooled via CO and SiO emission. In SN 1987A 
the mass of cool oxygen in the CO dominated region is 
estimated as 0.2 Mq (Liu & Dalgarno 1993). With a com- 
parable oxygen mass hidden in the SiO region we thus 
miss about 0.4 A4q of oxygen. Therefore the total oxygen 
mass must be « 1.6 Mq in rough agreement both with 
the estimate from the HS T spec trum at nonthermal exci- 
tation phase (Chugai et al. 1997) and predictions of stellar 



evolutio n mo dels for an 18-22 A Iq pro genitor (Woosley 



Weaver 1995; Thielemann et al. 1996) 



Omitting details, we conclude that the test of nebular 
model in the case of SN 1987A is successful and demon- 
strates that the model is able to recover reliable values of 
important parameters. 

3.3. SN 1997D: low and high-mass models 

Now we turn to the nebular spectrum of SN 1997D at 
t « 300 d. First, the 6 Mq case based on the hydrody- 
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Fig. 9. Computed (thick hne) and observed (thin hne) 
emission hnes in SN 1997D spectrum on day 296. Model 
Ml (a) with primordial metal abundance 0.3 of cosmic 
describes [Ca II] doublet better than the model M2 (b) 
with the cosmic abundance. 



namical model (Section 2.4) will be considered. Some re- 
finement of the hydrodynamical model is needed, however, 
to apply it to nebular epoch. The amount of metals in the 
mixed core is specified assuming that masses of metals and 
oxygen are equal. This is a reasonable assumption for a 
low-mass pre-SN. The oxygen abundance in He-rich mat- 
ter (component B) was assumed one tenth cosmic, while 
a carbon abundance of 0.03 is adopted for He-rich mat- 
ter in the 6 Mq model. We then also consider the 24 Mq 



case based on the model by Turatto et al. (1998) with the 



composition taken from Nomoto & Hashimoto (1988) 



The low-mass nebular model of SN 1997D fits the ob- 
served spectrum fairly well (Fig. m) with the optimal 
choice of parameters represented by model Ml (Table |l|). 
The [Ca II] 7300 A profile was reproduced for the core 
velocity Vc — 600 ± 30 km s^^. This parameter is of pri- 
mary importance, since it determines the absolute mass 
of the core components with the adopted density struc- 
ture. We failed to fit the absolute fiux of this line with a 
cosmic primordial abundance adopted for the model M2 
(Table H and Fig. 0b), while the primordial abundance 0.3 
of cosmic in model Ml provides an excellent fit. The oxy- 
gen doublet intensity is determined primarily by the mass 
of the 0-rich matter, although some 20% come from He 
and H-rich matter. The value of 0.035 Mq is corrected for 
the unseen cool oxygen assuming that we see 3/4 of all 
the pure oxygen in the [O I] doublet as in SN 1987A. A 
strong oxygen overdensity is not required. We found that 



Fig. 10. Computed nebular spectrum (thick line) in the 
high-mass case along with the observed (thin line) emis- 
sion lines in SN 1997D spectrum on day 296. (a) - model 
M3 without reduction of oxygen line-emitting mass; (b) - 
the same model but with five-fold reduction of the oxygen 
line-emitting mass. Neither of the model spectra fits the 
data. 



Xo = 1-2 in model Ml provides somewhat better agree- 
ment with the observed ratio of [O I] doublet components 
than Xo = 1- The amount of He- rich matter is a lower 
limit for the mass of the He shell in the pre-SN. One may 
admit up to 1 Mq of helium mixed microscopically with 
the H-rich envelope without a notable effect on line inten- 
sities. 

Unfortunately our nebular model is not applicable to 
the earlier nebular spectrum of SN 1997D on day 150. 
The reason is the significant optical depth in the Paschen 
continuum predicted by the model. In such a situation the 
multilevel statistical equilibrium must be solved together 
with a full radiation transfer, which is beyond the scope 
of this paper. Moreover, we found that the observed Ha 
profile at this epoch is odd exhibiting a significant redshift 
of unclear origin. A prima face explanation assuming ^^Ni 
asymmetry cannot accommodate to the late time nebular 
spectrum (« 300 d) lacking such an asymmetry. 

To evaluate the uncertainty related with the assump- 
tion of the same fraction (1/4) of cool 0-rich gas as in 
SN 1987A, we compared parameters relevant to molecular 
formation (density and temperature) in the 0-rich mat- 
ter at a similar nebular epoch. We found that the density 
of 0-rich gas in SN 1997D is lower, while the temper- 
ature is somewhat higher compared to SN 1987A. Both 
parameters suggest therefore that formation of molecules 
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estimated ^^Ni mass is about 0.002 M© in accordance with 



52 - 
51 - 

50 - 



-1 
-2 
-3 







a 


- 


93J 


87A 




- 


^^ 


^ 




- 


97D 






- 


93J 


87A 


b 


- 


^^ 


^^ 






97D 






- 


^^ 






■ 




1 


1 





10 



20 



30 



Fig. 11. Kinetic energy (a) and ^^Ni mass (b) as a func- 
tion of progenitor mass for three core-collapse supernovae. 

in SN 1997D cannot be more efficient than in SN 1987A, 
which means that the fraction of unseen pure oxygen in 
SN 1997D does not exceed that in SN 1987A. With the un- 
certainty of the core velocity the estimated range of pure 
oxygen mass in SN 1997D is 0.02-0.07 Mq. 

We applied the nebular model to the high-mass case 
(M = 2AMq). Due to the large mass of the He/0 core 
the velocity of the core boundary is too high and incon- 
sistent with the observed [Ca II] doublet profile. Mixing 
all the freshly synthesized helium with the hydrogen en- 
velope reduces the core velocity but not sufficiently to re- 
solve this controversy (Fig. [lO|a). Another serious problem 
is a high [O I] doublet flux and wrong doublet ratio. The 
five-fold reduction of amount of line-emitting oxygen, pre- 
sumably due to molecular formation and cooling, allevi- 
ates the problem of total flux in the [O I] doublet. Yet 
the problem of high /(6364)//(6300) ratio in this model 
remains (Fig. 10 3). 

Summing up, we found that the hydrodynamical model 
of moderate mass ejecta (6 Mq) which contains low 
amount of freshly synthesized oxygen (0.02-0.07 Mq) is 
consistent with the nebular spectra of SN 1997D. The 
model of high-mass ejecta as such is incompatible with 
the observed nebular spectra. 



4. Discussion 

We attributed SN 1997D to the SN IFF event character- 
ized by the kinetic energy E k, 10^° erg and ejecta mass 
6 ± \Mq. The ejecta are dominated by H-rich matter and 
contains 0.02-0.07 Mq of freshly synthesized oxygen. The 



the value found by Turatto et al. ( |1998D . The pre-SN had 
a moderate radius of 85 Rq and possibly a low primordial 
metallicity, 0.3 cosmic. 

At first glance the suggested low metallicity of SN 
1997D disagrees with the assumed cosmic abundance of 
barium which implies a relative overabundance by a factor 
of three. The latter should not confuse us, however, after 
SN 1987A in which the relative barium overabundance is 



about five (Mazzali et al. 1992) for the comparable metal- 
licity of both supernovae. Amazingly, the relatively small 
pre-SN radius and low primordial metallicity of SN 1997D 
both are reminiscent of SN 1987A. Fossibly it reflects some 
trend for low metallicity progenitors to have smaller pre- 
SN radii compared to SNe II-P with cosmic metallicity. 

When combined the ejecta mass and the collapsed 
core mass (presumably 1.4 Mq), the total pre-SN mass 
amounts to 6-9 M© prior to outburst. The main-sequence 
progenitor likely was more massive because of a possible 
wind mass-loss. In the context of general results of stellar 
evolution theory, the low mass of freshly synthesized oxy- 
gen (< 0.1 Mq) is compelling evidence that the progenitor 
of SN 1997D was a main-sequence star from the 8-12 Mq 
range. These stars are known to end their life with very 
low amount (< 0.1 Mq) of synthesized oxygen (Nomoto 
198^ Woosley |l986| ). Ejecta of CCSN produced by such 
stars must contain very small amount of ^^Ni, significantly 



less than normal CCSN (Woosley 198(:), which is also in 
line with SN 1997D. 

The fact that at least some CCSN originating from the 
8-12 Mq mass stars have low kinetic energy [E « 10^" 
erg) and eject small amounts of ^^Ni (« 0.002AfQ) mod- 
ifies a picture of CCSN with "standard" kinetic energy 
of - lO'^i erg and ^^^Ni mass of 0.07-0.1 Mq. The new 
situation in the systematics of CCSN is visualized by 
the E - Afnis and ^^Ni mass-Mms plots (Fig. |l^) which 
show the position of SN 1997D along with two other well 
studied CCSN, SN 198 7A (W oosley p^ ; Shygeyama & 
Nomo to |l990t Utrobin |1993D and SN 1993 J (Bar tunov 
et al. |1994|; S higeyama et al. |1994| ; Woosley et al. |1994 ; 
Utrobin 1996). The primary significance of this plot is that 
both nearly constant kinetic energy {^ 10°^^ erg) and ^^Ni 
mass (~ O.O8M0) in the range of progenitor masses be- 
tween a: 13M0 and « 2OM0 abruptly drop at the low end 
of massive star range producing CCSN (around 10 Mq). 
It would be not unreasonable to consider that SN 1997D 
is a prototype for a new family of CCSN (below referred 
to as "weak CCSN") which occupies the same place on 
the E-M^s and ^'^Ni mass-M^g plots as SN 1997D. 

Unfortunately, there are no clear theoretical predic- 
tions in regard to weak CCSN. Yet current trends in 
the core-collapse modelling seem to be generally consis- 
tent with Fig. 11, Two explosion mechanisms are related 
to producing CCSN: prompt (core rebound) and delayed 
(neutrino-driven mechanism). For 8-10 Mq progenitors 
the prompt mechanism attains its highest efficiency (Hille- 
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brandt et al. 1984) with the kinetic energy of ejecta ~ lO'^" 
erg (Baron & Cooperstein 199C), while the delayed mech- 
anism, on the contrary, has the lowest efficiency in this 
mass range with similar energy ~ 10^" erg (Wilson et 
al. 1986). Thus both explosion mechanisms remain viable 



terpart of a galactic, weak CCSN is SN 1181. With the 
absolute luminosity of —13.8 mag at maximum (Green & 



in the context of SN 1997D. However, possibly only the 
neutrino-driven mechanism is able to account for the ki- 
netic energy increase with the progenitor mass in the range 
from abou t 10 M© to « ISM© (Wilson et al. |1986| Bur- 
rows |l99l). 



How frequent are SN 1997D-like phenomena? The first 
thought is that they are extremely rare, since among ^ 10^ 
identified SN 11-P events only one such case has been dis- 
covered so far. However, with the low absolute luminosity 
(« —14 mag) and the brief plateau duration (40-50 days) 
compared to normal SN II-P characteristics (~ —16.5 mag 
and 80-100 days, respectively) it would not be surprising, 
if SN 1997D-like events were as frequent as ~ 20% of nor- 
mal SN II-P rate. Such a rate might be maintained by 
progenitors from the mass interval AMms « IMq in the 
vicinity of main-sequence mass w IOA/q. 

Progenitors from the 8-12 Mq mass range were sug- 
gested earlier as counterparts for supernovae with a dense 
circumstellar wind, low ejecta mass (~ ^Mq), and pos- 
sibly normal kinetic energy (e.g. SN 1988Z, Chugai & 
Danziger [1994D. The present attribution of SN 1997D to 



the same mass range introduces some dissonance with the 
former conjecture. In reality, this controversy is not serious 
since 8-12 Mq progenitors are characterized by very com- 



plicate evolutions of their cores (Nomoto 1984; Woosley 



1986) and therefore a different outcome for slightly dif- 



ferent initial mass is quite conceivable. Moreover, it may 
well be that with a normal metallicity presupernova of 
weak CCSN also vigorously loses mass and explodes in a 
dense wind thus producing a luminous supernova (possi- 
bly SN Iln) due to the ejecta wind interaction. 

Another intriguing possibility is that a presupernova 
of weak CCSN might lose all the hydrogen envelope in a 
close binary system before the explosion. In this case weak 
CCSN will be a mini-version of SN lb with low explosion 
energy, low amount of ^^Ni, and, eventually, low luminos- 
ity. Unfortunately it will not be easy to detect such events. 

If weak CCSN are as frequent as ~ 20% of all CCSN, 
a good fraction of galactic population of SNR may be 
related to these supernovae. We cannot miss the oppor- 
tunity to speculate that at least two historical super- 
novae SN 1054 and SN 1181 might be identified with weak 



CCSN. Nomoto ( |1984D already argued that the Crab Neb- 
ula was created by CCSN with the progenitor mass around 
9 Mq. The luminosity of SN 1054, which was normal for 
SN 11, could be explained then by the interaction of ejecta 
with a dense pre-SN wind. This suggestion in fact is a 
modification of the earlier idea that the initial phase of 
the light curve of SN 1054 could be powered by the shock 
wave propagating in the circum stellar (r r^ 10"'^^ cm) enve- 
lope (Weaver & Woosley 1979 ) . The second possible coun- 



GuU |1982| ) and half year period of visibility SN 1181 is 
the closest analogue of SN 1997D. Low radial velocities 
(< 1100 km s-i) of SN 1 181 filaments claimed by Fesen, 
Kirshner & Becker (1988) seem to strengthen this identi- 
fication. If the association of SN 1181 is correct then we 
expect to find very low amount of newly synthesized oxy- 
gen and iron-peak elements in this supernova remnant. 
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